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The dynamical mean-field theory (DMFT) combined with the fluctuation exchange (FLEX) 
method, namely FLEX+DMFT, is an approach for correlated electron systems to incorporate both 
local and nonlocal long-range correlations in a self-consistent manner. We formulate FLEX+DMFT 
in a systematic way starting from a Luttinger-Ward functional, and apply it to study the d- wave 
superconductivity in the two-dimensional repulsive Hubbard model. The critical temperature (T c ) 
curve obtained in the FLEX+DMFT exhibits a dome structure as a function of the filling, which 
has not been clearly observed in the FLEX approach alone. We trace back the origin of the dome to 
the local vertex correction from DMFT that renders a filling dependence in the FLEX self-energy. 

We compare the results with those of GW+DMFT, where the T c -dome structure is qualitatively 
reproduced due to the same vertex correction effect, but a crucial difference from FLEX+DMFT is 
that T c is always estimated below the Neel temperature in GW+DMFT. The single-particle spec¬ 
tral function obtained with FLEX+DMFT exhibits a double-peak structure as a precursor of the 
Hubbard bands at temperatures above T c . 

PACS numbers: 71.10.Fd, 74.20.-z, 74.25.Dw 


I. INTRODUCTION 

Despite a long history of physics of the high-T c 
cuprate, 1,2 we are still some way from a full understand¬ 
ing of the superconductivity. There is a general con¬ 
sensus that the supercurrent flows on each Cu-0 plane, 
which can be modeled by the repulsive Hubbard model 
on the square lattice. There are actually two essen¬ 
tial factors here: the repulsive Hubbard interaction can 
give rise to a pairing interaction in the d -wave channel 
mediated by antiferromagnetic spin fluctuations, 3 while 
the very same interaction also introduces Mott’s metal- 
insulator transition 4 that hinders the superconductivity 
around half-filling for strong enough interactions. Cap¬ 
turing these two features simultaneously still remains a 
theoretically challenging task. As numerical methods for 
treating the strongly correlated electron systems, there 
are the exact diagonalization and quantum Monte Carlo 
(QMC) methods 5 that are exact within numerical errors, 
but the former can only deal with limited system sizes, 
while the latter suffers from the sign problem. 

However, we do have theoretical methods that can deal 
with each of the d-wave pairing and Mott’s transition sep¬ 
arately: Namely, we have on one hand the fluctuation- 
exchange (FLEX) approximation, 6 one of the perturba¬ 
tive methods for many-body physics that can describe 
the spin-fluctuation mediated d-wave pairing. On the 
other hand, we have the dynamical mean-field theory 
(DMFT), 7-9 which can describe the Mott transition. To 
be more precise, the FLEX describes the momentum de¬ 
pendence of the effective pairing interaction mediated by 
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the antiferromagnetic spin fluctuations, which is essen¬ 
tial for the anisotropic pairing, 6 but the method, being 
perturbative, cannot describe the Mott transition in the 
regime close to the half-filling. The DMFT, although 
mean-field theoretic, describes Mott’s insulator in terms 
of the (non-perturbative) correlation effect that is local 
(i.e., momentum-independent) but dynamical (i.e., incor¬ 
porating temporal fluctuations), and becomes exact in 
the limit of infinite spatial dimensions of a lattice model. 7 

There are many extensions of DMFT to include mo¬ 
mentum dependence of the self-energy. 10-24 One is the 
cluster extension of DMFT, 10,11 which is employed, 
e.g., for explaining the pseudogap in the cuprates as 
a momentum-selective Mott transition. 12,13 However, in 
practice it is quite hard in this scheme to attain large clus¬ 
ter sizes and to incorporate spatially long-ranged compo¬ 
nents in the self-energy in a strongly correlated regime. 
It is also computationally very demanding to treat the d- 
wave superconducting phase, or to extend to more com¬ 
plicated models such as multi-orbital systems with a large 
cluster size retained in the cluster DMFT. 

More realistically, we have alternative and numeri¬ 
cally feasible extensions of DMFT that combines DMFT 
with a certain resummation technique of nonlocal self¬ 
energy diagrams, such as GW+(E)DMFT, 14,15 DTA, 16,17 
and the dual-fermion approach. 18 20 These schemes can 
treat momentum-dependent self-energies describing non¬ 
local long-range correlations with some selected diagrams 
taken into account. This has motivated us to take the 
FLEX+DMFT method, 24,25 where nonlocal FLEX dia¬ 
grams are considered on top of DMFT local diagrams for 
the self-energy. We have opted for a method that evokes 
FLEX among other diagrammatic methods, since we are 
interested in the d-wave superconductivity mediated by 
antiferromagnetic fluctuations, which can be explicitly 
treated with FLEX. 
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In the present paper, we extend the FLEX+DMFT 
method to deal with the d-wave superconductivity in 
the two-dimensional repulsive Hubbard model, while a 
FLEX+DMFT has been applied to the normal phases 
of the Hubbard model in Ref. 24. To this end, we con¬ 
struct the Luttinger-Ward functional for FLEX+DMFT, 
where double counting of local diagrams from FLEX and 
DMFT parts is unambiguously subtracted. Starting from 
the Luttinger-Ward functional formalism guarantees the 
conserved nature of DMFT (as well as FLEX) retained, 
which is not always the case with other diagrammatic 
extensions of DMFT. We then apply this FLEX+DMFT 
to the d-wave superconductivity in the two-dimensional 
repulsive Hubbard model to obtain the superconducting 
phase diagram. 

We find that the FLEX+DMFT result exhibits a T c - 
dome structure of the superconducting phase diagram 
against band filling, which has not been observed in 
FLEX alone. We identify the origin of the dome to 
the local vertex correction from DMFT that renders a 
filling dependence in the FLEX self-energy. To elabo¬ 
rate this point, we compare this with the GW+DMFT 
method, in which only bubble diagrams are used to ex¬ 
tend DMFT in considering a nonlocal self-energy correc¬ 
tion, whereas both bubbles and ladders are included in 
FLEX+DMFT. The GW+DMFT result also exhibits a 
T c -dome structure, but, unlike the FLEX+DMFT result, 
T c in GW+DMFT is always below the Neel temperature, 
i.e., the antiferromagnetic order dominates over d-wave 
superconductivity for the whole filling range. We have 
also obtained the single-particle spectral function with 
the FLEX+DMFT, which exhibits a double-peak struc¬ 
ture above T c with a precursor of the Hubbard bands. 

While the present scheme does not consider vertex 
corrections to the nonlocal ladder diagrams unlike the 
dual-fermion approach which is recently applied 26 to the 
superconductivity in the Hubbard model, an advantage 
of the present method is that we define the Luttinger- 
Ward framework, which enables us to treat the normal 
self-energy and the anomalous (d-wave) self-energy on an 
equal footing as derivatives of the same Luttinger-Ward 
functional. 


II. FLEX+DMFT FUNCTIONAL 

Let us formulate the FLEX+DMFT method by in¬ 
troducing a Luttinger-Ward functional <f>, 27 which ba¬ 
sically consists of FLEX and DMFT diagrams. How¬ 
ever, there is a double counting of local self-energy di¬ 
agrams between the two contributions, which must be 
subtracted. We show that the double counting term is 
uniquely identified if one demands the conserving nature 
of the formalism. Namely, we regard each of DMFT 
and FLEX as an approximation for the exact Luttinger- 
Ward functional of the dressed Green’s function G to 
propose a new functional, in a manner similar to the 
GW+(E)DMFT scheme. 14,15 In DMFT, the approximate 


functional, 4>dmft, is the sum of all types of the ring di¬ 
agrams that only contain the local Green’s function G\ oc . 
On the other hand, the approximate functional in FLEX, 
^flex, is the sum of specific (bubble and ladder) dia¬ 
grams as shown in Fig. 1(a), which basically correspond 
to spin and charge fluctuations. 

Then we can propose a functional in the FLEX+ 
DMFT scheme as 


^FLEX+DMFT [G] = ^DMFt[Gi oc ]+^FLEx[G]—$flEx[^1oc] } 

(1) 

where we have subtracted the local part of the FLEX 
functional $flex[ g i°c] with Gi oc = (l/Vt) J2k G ( k ) 
(TVfc: number of k points) to avoid the double counting. 
Since both 4>dmft and 4>flex are expressed as function¬ 
als of dressed Green’s functions, the overlap between the 
two is uniquely determined as a set of diagrams in the 
^flex[G] that only contain local dressed Green’s func¬ 
tions, which is nothing but 4 >flex[£uqc]- We then obtain 
the self-energy in this scheme as a functional derivative, 


^flex+dmft[G] = 
53>DMFt[Gi, 


£Gi 0 


+ 


(^FLEX+DMFT 

SG 

<5($FLEX - $FLEx) 

SG 


( 2 ) 


This way we retain the conserving nature of the 
approximation. 28,29 The first term in the last line 
of Eq. (2) is the local DMFT self-energy Fi mp = 
£$dmft[Gioc]/£Gio C . The second term in Eq. (2), 

V'nonloc _ <K^FLEX “ ^FLEx) 

^FLEX = Jq 

= S FL ex[G] - SpLExt^loc], (3) 


is the difference between the FLEX self-energy con¬ 
structed from the lattice Green’s function G and that 
from the local Green’s function G\ oc . Note that Fflex 
contains some contributions from local parts of the self¬ 
energy, i.e., FpLEx+ 7^ 0 (i: label of lattice sites). For 
example, SpLEXii contains a diagram displayed in the 
left-hand side of Fig. 1(b) with i ^ j, while a diagram 
shown on the right-hand side of Fig. 1(b) does not belong 

v non l° c 
ro ^flex+- 

The self-consistency loop, which has to be a double 
loop in the present combined scheme, is depicted in 
Fig. 1(c): To start with, we define the DMFT mapping of 
a lattice model to an impurity model in such a way that 
Green’s function of the mapped impurity model, Gi mp , 
coincides with the local Green’s function for the original 
lattice model, G\ oc . The local self-energy Fi mp is calcu¬ 
lated in the DMFT part of the self-consistency loop. The 
nonlocal part of the self-energy Xflkx 1S ^ en calculated 
in the FLEX loop. We combine both Fi mp and Splex 
to obtain the full self-energy, from which we construct 
new (full and local) Green’s functions. We update each 
of them (E imp , E^ex) alternately by using the corre¬ 
sponding loops until the whole loops [Fig. 1(c)] converge. 
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FIG. 1: (Color online) (a) The Luttinger-Ward functional for 
FLEX. 30 (b) An example of a diagram contained in Xflex,^ 
(left, with i ^ j) and omitted (right), (c) Self-consistent 
double loops in the FLEX+DMFT formalism. The algorithm 
flow for the loops is indicated by the red arrow. 


The present scheme may be viewed as a new diagram¬ 
matic extension of the DMFT that incorporates vertex 
corrections into the (local part of) FLEX scheme. FLEX 
itself, being a perturbative method, is considered to be¬ 
come exact in the weak-coupling limit, while DMFT be¬ 
comes exact in the atomic limit. Since FLEX+DMFT 
formalism here incorporates the functionals that domi¬ 
nate in either limit, it is expected to describe spin fluc¬ 
tuation effects and Mott’s physics simultaneously. 


III. APPLICATION TO THE 2D HUBBARD 
MODEL 


Let us apply the FLEX+DMFT method to the repul¬ 
sive Hubbard model on the square lattice, with a Hamil¬ 
tonian, 


U = Y1 e ( k ) C l,a C k,a + U Y1 n *,t n U- 

k,cr i 


Here a creates an electron in a Bloch state with wave- 
vector k = ( k x , k y ) and spin cr, U is the on-site repulsion, 
and = c\ a Ci j(T is the number operator. The two- 
dimensional band dispersion is given as 


e(fe) = —2£(cos£+ + cos k y ) 

— At'cosk x cosk y — 2t"(cos2k x + cos2 k y ) — /q (5) 

where £, P, and t" represent the nearest-neighbor, second- 
neighbor, and third-neighbor hoppings, respectively, 
while /i is the chemical potential. We shall compare 
the case with the nearest-neighbor hopping only (t'/t = 
t"/t = 0) with the case of t'/t = —0.20, t"/t = 0.16, 
which are the values estimated for a typical hole-doped, 
single-layered cuprate, HgBa 2 Cu 04 + <5 with T c ~ 90K, 
with first-principles methods. 31,32 Hereafter we take \t\ 
as the unit of energy. 

In the single-band Hubbard model, the FLEX self¬ 
energy is computed as 


^FLEx(fc) = 


1 y^r 3 7- r 2 Xo{k - k') 

N k p^V2 1 — U xo(fc — k') 


, 1 U 2 Xo (k-k') 

2 1 + Uxo(k — k') 


-U 2 Xo (k-k') G(k'), 


(6) 


where f3 is the inverse temperature, k = (o; n , k) with uo n 
the fermionic Matsubara frequency, G(k) is the Green’s 
function, and 


Xo(q) ^~r2-p^2G(k + q)G(k) (7) 

k 


is the irreducible susceptibility. We can calculate Eflex 
by replacing G with G\ oc in Eqs. (6) and (7). 

To obtain Xi mp in the DMFT procedure, we have 
to solve the impurity problem in DMFT. Among var¬ 
ious impurity solvers, here we adopt the modified it¬ 
erative perturbation theory (modified IPT), where the 
original IPT is modified for systems without particle- 
hole symmetry. 33 The method is not computationally de¬ 
manding, which facilitates a scanning over a wide param¬ 
eter region to obtain the phase diagram, and also enables 
us to approach a region with large antiferromagnetic fluc¬ 
tuations where FLEX convergence critically slows down. 
We have confirmed for various values of parameters 
that the continuous-time quantum Monte Carlo impurity 
solver 34,35 implemented with the ALPS library 36,37 gives 
similar values for the eigenvalue of Eliashberg’s equation 
even away from the half-filling. 

When Green’s function is obtained, we plug it into the 
linearized Eliashberg equation, 

XA(k) = -^ X V *( k ~ k')G{k')G{-k')A{k'). (8) 


Here A (k) is the anomalous self-energy, which is the gap 
function up to the renormalization factor, and 


V eS (k) = u + ^u 2 


Xo(fc) 

1 - U X o(k) 


+2 Xo(fc) 

2 1 + U X o(k) 


(4) 


(9) 
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FIG. 2: Effective pairing interaction in the FLEX. 


is the effective pairing interaction (Fig. 2), where A is 
the eigenvalue for Eliashberg’s equation, with supercon¬ 
ducting transition identified as the temperature at which 
A = 1. 

At this point we should mention about the consis¬ 
tency of the approximate functional form $flex+dmft 
and the linearized Eliashberg equation, Eq. (8). The 
Luttinger-Ward functional can be extended to incor¬ 
porate the anomalous part, and the extended func¬ 
tional <F[G, F^F] is related to the anomalous self-energy 
through A = 8$/8F\ where F ^ is the anomalous Green’s 
function, for which we should consider the local correc¬ 
tion to the anomalous self-energy A as Aflex+dmft = 
Aflex + Aioc as in Eq. (2) for the normal self-energy. 
Now, our interest here is the anisotropic, d-wave pair¬ 
ing instability in the repulsive model, for which we can 
ignore the local correction to the anomalous self-energy 
Aioc which does not depend on momentum. The remain¬ 
ing term Aflex = £4>flex[G, F\ F]/8F^ is the same as 
the right-hand side of the linearized Eliashberg equation 
(8) if we linearize the anomalous part. 38 Then our formal¬ 
ism treats the normal and anomalous self-energies con¬ 
sistently, as functional derivatives of the same Luttinger- 
Ward functional 4>flex+dmft- This is an advantage of 
using the Luttinger-Ward functional formalism in con¬ 
structing a new scheme. 




IV. RESULTS 

We show the superconducting phase diagram of 
the two-dimensional Hubbard model obtained in the 
FLEX+DMFT in Fig. 3, right panels, where the FLEX 
result is also displayed in the left panels for comparison. 
We can immediately see that T c exhibits a dome structure 
in the FLEX+DMFT. This sharply contrasts with the 
FLEX result, where T c has been known to almost mono- 
tonically increase toward half-filling with some rounding 
off 39 The p resence Q f the T c dome in the FLEX+DMFT 
and its absence in the FLEX are seen for both the simple 
square lattice with t' = t" = 0 [Fig. 3(a)] and the case of 
t' = —0.20, t" = 0.16 [Fig. 3(b)]. 


FIG. 3: (Color online) Phase diagram against the band fill¬ 
ing n and the temperature T/t in the FLEX+DMFT (right 
panels) as compared with FLEX (left). Here we take U/t = 
4.0, 5.0 and (a) t'/t — t"/t = 0, or (b) t'/t = —0.20, t"/t = 
0.16. We also plot the AF phase boundaries (dotted lines) in 
the normal region, while color shading highlights the super¬ 
conducting region with T c > Taf at each filling. 


For the simple square lattice (t f = t" = 0), we cannot 
approach a region very close to half-filling because the 
antiferromagnetic (AF) fluctuations prevent the FLEX 
self-consistency loop from converging. For the same rea¬ 
son, it is difficult to attain convergence for systems with 
larger U. As a measure of the AF order, we evaluated 
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FIG. 4: (Color online) (a) Comparison between the FLEX 
local self-energy ImS^Ex (red plus signs) and the DMFT 
self-energy ImXi mp (green crosses) for the filling n = 0.70 
(underdoped; left panel), 0.88 (optimally doped; center), and 
1.0 (half-filled; right), (b) The difference, ImXpL EX — ImXi mp , 
for n — 0.70 (red crosses), 0.88 (green circles), and 1.0 (blue 
squares). Here we take U/t = 4.0, ft = 20, t'/t = —0.20, and 
t"/t = 0.16. 


the AF phase boundaries (dashed lines in Fig. 3) de¬ 
termined from maxk[Uxo(k)](= 0.99 here), 40 which is 
usually adopted in FLEX-type schemes to take account 
of the effect of the quasi-two-dimensional nature (e.g., 
in three-dimensional layered systems), although FLEX- 
type approaches are known to obey the Mermin-Wagner 
theorem that forbids finite-temperature AF phase tran¬ 
sitions in an isolated two-dimensional system. The esti¬ 
mated AF transition temperature Taf becomes higher 
than the superconducting T c as one approaches half¬ 
filling as shown in Fig. 3, where the color shaded region 
indicates the superconducting phase with T c > Taf (be., 
superconductivity dominating antiferromagnetism). The 
result suggests that a part of the T c dome is taken over by 
the AF phase in the case of t' = t" = 0 [Fig. 3(a), right] 
and t' = —0.20,t" = 0.16, V — 5 [Fig. 3(b), right]. For 
a smaller U = 4, by contrast, we have an almost full T c 
dome with T c > Taf for t' = —0.20 ,t" = 0.16 [Fig. 3(b), 
right]. These are a key result in the present work. 

Now let us identify the physical origin of the ap¬ 
pearance of the T c dome in the FLEX+DMFT. In 
FLEX+DMFT, the self-energy is obtained from the 
FLEX and DMFT self-energies as Xflex+dmft = 

SfLEX - (SpLEX “ Sim p) [ EqS - and ( 3 )]’ Le -’ a P art 


of the local self-energy is replaced from that in FLEX 
with that in DMFT. Thus the quantity X eeex — ^im P 
represents the difference in the self-energy effect be¬ 
tween FLEX and FLEX+DMFT. We can actually take 
a look at X eeex and Xi mp , with fillings n = 0.7 (un¬ 
derdoped), 0.88 (optimally doped), and 1.0 (half-filled) 
in Fig. 4 [the parameters are taken to be U = 4.0, /? = 
20, t' = —0.20 ,t" = 0.16, which corresponds to Fig. 3(b), 
right panel]. We first notice that the magnitude of the 
DMFT self-energy X imp is smaller than that of FLEX 
^flex’ which means that the overestimation of the self¬ 
energy generally known to exist in FLEX is remedied 
in FLEX+DMFT by the DMFT (local) vertex correc¬ 
tions. More importantly, we can see that the differ¬ 
ence, ImXpp EX — ImX imp [Fig. 4(b)], has a clear fill¬ 
ing dependence, and it increases with doping. Since 
^FLEX+DMFT = ^FLEX ~ (^FLEX — ^imp), the result 
indicates that the reduction of the FLEX+DMFT self¬ 
energy due to DMFT correction is reduced as one ap¬ 
proaches half-filling. Thus T c tends to be suppressed near 
half-filling as compared to that of FLEX because of the 
filling-dependent self-energy reduction in FLEX+DMFT. 
On the other hand, the pairing interaction itself aris¬ 
ing from spin fluctuations becomes stronger toward half¬ 
filling due to better band nesting, as reflected in the 
FLEX result [Fig. 3, left panels] with T c almost mono- 
tonically increasing toward half-filling. Therefore, the 
FLEX+DMFT contains two factors with opposite fill¬ 
ing dependencies, and we conclude that the T c dome 
in FLEX+DMFT arises from the combined effect of the 
nesting and filling-dependent self-energy reduction. 

In the FLEX+DMFT scheme, the self-energy reduc¬ 
tion from DMFT takes place only in the local part, 
while the nonlocal self-energy is still considered to be 
overestimated, especially for ladder diagrams. 24 To ex¬ 
amine this effect, we compare the present method with 
GW+DMFT, where only the bubble diagrams are con¬ 
sidered for the self-energy and the pairing interaction 
(whereas both bubbles and ladders are included in FLEX 
and FLEX+DMFT). We show the GW+DMFT phase 
diagram, along with the GW result for comparison, in 
Fig. 5 for t' = —0.20 and t n = 0.16. We can see 
that, although the T c dome structure remains in the 
GW+DMFT result, T c is much reduced from the result 
of FLEX+DMFT. On the other hand, the AF transition 
temperature is much higher in GW+DMFT than that of 
FLEX+DMFT. This makes the region in the dome where 
T c > Taf [highlighted with color shadings in Fig. 5(a)] 
very narrow in GW+DMFT. In fact, for t' = t" = 0 the 
AF instability becomes so strong that we cannot even 
obtain superconducting phase boundaries for the whole 
region of the fillings considered. 

In Fig. 5(b), we display the GW local self-energy 
ImX 1 ^ as compared with the DMFT self-energy IrnXi mp 
for the filling n = 0.70,0.88,1.0 with U = 4.0 and 
P = 50. We can see that the filling dependence is simi¬ 
lar to those in the FLEX+DMFT in that the difference 
between the two self-energies increases with the doping. 
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FIG. 5: (Color online) (a) Phase diagram against temperature 
T/t and band filling n in GW+DMFT approximation (right) 
as compared with that in GW approximation (left) for U/t — 
4.0 (green squares) or U/t = 5.0 (red circles), (b) Comparison 
between the GW local self-energy IiuEgw (red plus signs) and 
the DMFT self-energy ImSi mp (green crosses) for the filling 
n = 0.70 (left panel), 0.88 (center), and 1.0 (right) with U/t = 
4.0 and (3 — 50. Here we take t'/t = —0.20, t"/t — 0.16. We 
also plot antiferromagnetic (AF) phase boundaries (dotted 
lines) in the normal region, while color shading highlights the 
dome in the region where T c is above the AF boundary. 


Hence we can conclude that the existence of the T c dome 
is not an artifact in FLEX+DMFT, but is robust in both 
FLEX+DMFT and GW+DMFT arising due to the same 
local vertex correction effect. The overestimation of non¬ 
local self-energy thus does not affect the existence of the 
T c dome itself. 

The reason that the magnitude of T c is much smaller 
in GW+DMFT than in FLEX+DMFT is because lad¬ 
der diagrams describing spin fluctuations are not taken 
into account in GW+DMFT. In this sense, GW+DMFT 
is closer to the mean-field theory than FLEX+DMFT, 



FIG. 6: (Color online) Double occupancy against band filling 
n in the FLEX+DMFT (red circles), DMFT (green crosses), 
FLEX (blue squares), GW+DMFT (light blue diamonds), 
and GW (purple triangles). Here we take U/t = 4.0, /3 — 20, 
t'/t = —0.20, and t"/t = 0.16. 


which is also reflected in the higher AF transition tem¬ 
perature in GW+DMFT. Concomitantly, the pairing in¬ 
teraction mediated by spin fluctuations is reduced, which 
acts to reduce the superconducting T c in GW+DMFT 
rather than in FLEX+DMFT. The fact that T c is always 
estimated below the Neel temperature in GW+DMFT 
suggests that GW+DMFT underestimates the spin fluc¬ 
tuation effect, and is not enough to describe the d-wave 
superconductivity mediated by spin fluctuations. Since 
the overestimated nonlocal self-energy in FLEX+DMFT 
is remedied in GW+DMFT, the accurate estimation of 
the spin-fluctuation effect is expected to lie between 
GW+DMFT and FLEX+DMFT. 

To see the strength of local correlation, we measure the 
double occupancy, 


++ = ^TY (EG). (10) 

We can see in Fig. 6 that the double occupancy becomes 
negative in the overdoped regime in FLEX, while this 
unphysical behavior is improved in FLEX+DMFT. We 
can regard this as another of the self-energy reduction 
effects: FLEX overestimates the correlation effect, while 
this is corrected by combining it with the DMFT. A simi¬ 
lar tendency is observed between GW and GW+DMFT, 
but the difference in the double occupancy is smaller. 
This should be because the self-energy reduction effect is 
smaller [see Figs. 4(a) and 5(b)]. 

Let us finally examine the spectral function, which 
is calculated via analytical continuation with the Pade 
approximation. The results for FLEX+DMFT, FLEX, 
GW, GW+DMFT, and DMFT at various fillings are 
shown in Fig. 7. We can see that the filling dependence is 
similar among FLEX, GW, and DMFT in that we have 
a single peak that slightly shifts and broadens as we ap¬ 
proach the half-filling. By contrast, in the method that 
combines DMFT with either FLEX or GW, the spectral 
function acquires a stronger filling dependence, where a 
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FIG. 7: (Color online) Spectral functions at filling n = 0.70 
(underdoped; dotted lines), 0.88 (optimally doped; green), 
1.0 (half-filling; red) in the FLEX+DMFT (top left) and 
GW+DMFT (bottom left) are compared with those in DMFT 
(top right), FLEX (middle right), and GW (bottom right). 
Here we take U/t = 4.0, f3 = 20, t'/t = —0.20, and 
t"/t = 0.16. 


marked double peak is observed at half-filling. Similar 
double-peak structures have been reported in the dual¬ 
fermion method 21 as an antiferromagnetic pseudogap, 
where the appearance of the double peak is consistent 
with the QMC result. Thus we can see that the interplay 
of the local and nonlocal long-range correlation effects in 
FLEX+DMFT and GW+DMFT gives the double peak, 
which is considered to be a precursor of the Hubbard 
bands with two peaks separated by about [/, while the 
system is metallic. 

If we look more closely at the momentum-resolved 
spectral function A(fe,u;), we observe that there is a re¬ 
gion in k space near the Fermi energy where the spec¬ 
tral weight becomes slightly negative. This might not 
be specific to the present method, since many exten¬ 
sions of DMFT do not guarantee positive-definite spec¬ 
tral weights. 20,22 Since the magnitude of the negative 
part is negligibly small (< 1%) in the present case and 
this tends to occur in the overdoped regime, this does 
not affect the phase diagram and the density of states 
[Figs. 3 and 7] in the underdoped regime. 

V. SUMMARY AND DISCUSSIONS 

We have employed the FLEX+DMFT approach in 
terms of a new Luttinger-Ward functional to study the 


superconductivity in correlated electron systems. This 
scheme is a diagrammatic extension of the DMFT, so 
that it can describe the d-wave superconductivity aris¬ 
ing from ^-dependent pairing interaction. The scheme, 
being formulated in terms of the Luttinger-Ward func¬ 
tional, also has a virtue of the normal and anomalous 
self-energies treated on an equal footing. We have ap¬ 
plied the FLEX+DMFT to the repulsive Hubbard model 
on the square lattice. We have found that FLEX+DMFT 
describes a T c dome structure, whose physical origin is 
traced back to a combination of opposite effects: The 
self-energy effect introduced in the FLEX+DMFT sup¬ 
presses the superconductivity more strongly toward the 
half-filling due to the local-correlation effect, while spin 
fluctuations become stronger toward the half-filling due 
to band nesting. We also compare the FLEX+DMFT re¬ 
sult with the GW+DMFT result, which reproduces the 
dome structure. This indicates that the dome is not 
an artifact of the overestimated nonlocal self-energy in 
FLEX+DMFT. 

Another observation is that there is a Pomeranchuk 
instability into electronic states with broken tetragonal 
symmetry in the case of t' = —0.20,^" = 0.16 in both 
FLEX+DMFT and GW+DMFT. In this case, solutions 
with the four-fold rotationally symmetric Fermi surface 
become unstable, where we end up with a solution that 
breaks this symmetry when we start the calculation from 
an asymmetric initial input. While this instability is in¬ 
teresting in its own right, we have concentrated on the 
symmetric case in this study, and leave the analysis of 
the Pomeranchuk instability to another publication. 

In order to improve the scheme to suppress the over¬ 
estimation of the nonlocal FLEX self-energy, we should 
consider the screening effect in the FLEX self-energy. For 
example, the two-particle self-consistent method 41 takes 
account of vertex correction effects by considering the 
sum rule for the susceptibility, while a similar technique 
is also used to reduce the overestimated spin fluctuations 
and their effect on the self-energy in DTA. 17 We expect 
these techniques bring some improvement to the present 
theory. 
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